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Abstract 

For two points p and q in the plane, a straight line h, called a highway, and a real v > 1, we define 
the travel time (also known as the city distance) from p and q to be the time needed to traverse a 
quickest path from p to q, where the distance is measured with speed v on h and with speed 1 in the 
underlying metric elsewhere. 

Given a set S of n points in the plane and a highway speed v, we consider the problem of finding 
a highway that minimizes the maximum travel time over all pairs of points in S. If the orientation of 
the highway is fixed, the optimal highway can be computed in linear time, both for the L\- and the 
Euclidean metric as the underlying metric. If arbitrary orientations are allowed, then the optimal 
highway can be computed in 0(n 2 log n) time. We also consider the problem of computing an optimal 
pair of highways, one being horizontal, one vertical. 

Keywords: geometric facility location, min-max-min problem, city metric, time metric, optimal highways 



1 Introduction 

Many facility location problems are min-max problems, where the task is to place a facility such that 
the maximum cost incurred by any customer is minimized. For example, if the customers are points in 
the plane and their cost of using a facility is the Euclidean distance to the facility, then the center of 
the smallest enclosing circle is the optimal facility location that minimizes the maximum distance to the 
customers. Recently, Cardinal and Langerman [7] introduced the class of min-max-min facility location 
problems, where customers have the choice of either using or not using the facility, and their cost is 
the minimum cost of these two options. Transport facility location problems are typical min-max-min 
problems: a new railway line or highway will not be used by a customer if using it does not improve the 
travel time compared to existing means of transportation. Cardinal and Langerman consider three such 
min-max-min problems, among them the following: Given a set P of pairs of points, find the highway (a 
straight line) that minimizes the maximum travel time over all pairs in P. Here, the travel time of a pair 
of points (a, b) is the time taken to travel from a to b, assuming constant speed anywhere in the plane, 
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Table 1: Overview of our results for one highway. 
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(1 + e)-approx. 


0(n 4 a(n)) 

0(n logn) 
0(n log(l/e)a(n) logn) 
0(n 2 Iog(l/e)a(n) logn) 



Table 2: Overview of our results for the highway cross (under the Li-metric). 

infinite speed along the highway, and travel to and from the highway parallel to the y-axis. Cardinal and 
Langerman show how to compute the optimal highway in expected time linear in the number of pairs. 

Previous work considering highways (also called transportation networks or roads) have focused on 
how to compute quickest paths among the customers (points) or their Voronoi diagrams under the metric 
induced by given highways. Abcllanas et al. [T] started work in this area, and discussed the Voronoi 
diagram of a point set given a horizontal highway under the Li-metric. They later also considered the 
problem under the Euclidean metric and studied shortest paths [2]. Aichholzcr ct al. [3] introduced the 
city metric induced by the Li-metric and a highway network that consists of a number of axis-parallel 
line segments. They gave an efficient algorithm for constructing the Voronoi diagram and a quickest- 
path map for a set of points given the city metric. The running time of their algorithms was recently 
improved by Gorke et al. [TT] and Bae et al. [BJ. Bae et al. [5] presented algorithms that compute the 
Voronoi diagram and shortest paths, using the Euclidean metric and more general highway networks 
whose segments can have arbitrary orientation and speed. They recently extended their approach to 
more general metrics including asymmetric convex distance functions [4] . 

In this paper, we, like Cardinal and Langerman, consider the problem of finding an optimal highway 
(a straight line) for a given set S of n points in the plane. We wish to place a highway such that 
the maximum travel time over all pairs of points (that is, the travel-time diameter of the point set) 
is minimized. Since there is a quadratic number of pairs, Cardinal and Langerman's algorithm takes 
expected quadratic time in this case. We show how to make use of the coherence between the pairs of 
points to get deterministic ncar-lincar-timc algorithms. 

We assume that we can travel anywhere in the plane with speed 1, and that we can travel along the 
highway with a given speed v > 1. We consider various versions of the problem. First, we distinguish 
whether the highway speed v is infinite (the easiest case) or finite. Second, we consider both the L\- 
metric (as in the city Voronoi diagram) and the L2-metric. In all cases, we show how to find the optimal 
highway with a given orientation (that is, the optimal horizontal highway). In the case of the Euclidean 
metric, we also consider how to find the optimal highway if we arc free to choose the orientation. (We 
note that chosing the orientation does not make sense for the Manhattan metric). Tabic [T] summarizes 
our results for all variations of the problem. 

We then consider the problem of placing a highway cross, that is, a pair of a horizontal and vertical 
highway. We can determine the optimal axis-aligned highway cross with infinite speed in O(nlogn) time, 
see Section [5] For constant speed the problem becomes considerably harder, even under the Li-metric. 
We give an exact 0(n 4 a(n))-time algorithm based on computing minima of upper envelopes. We also 
consider approximative solutions — all our results are summarized in Tabic [5] 

Throughout the paper we assume that the input point set S contains at least three points (if \S\ < 3, 
it is trivial to find an optimal highway) . 
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(a) horizontal (b) arbitrary orientation 

Fig. 1: Optimal infinite-speed highways (solid lines) and corresponding enclosing figures (shaded). 



Fig. 2: Optimal constant-speed horizontal highway (solid lines) and corresponding enclosing figure (shaded). 

2 The optimal highway for infinite speed 

As a warm-up exercise, let us consider the problem of finding the optimal placement of a horizontal 
highway, assuming the highway speed is infinite. 

Theorem 1 Given n points in the plane, the middle line of the smallest enclosing horizontal strip is an 
optimal horizontal highway of infinite speed. It can be computed in linear time. 

The easy proof is left to the reader. What is interesting is that the optimal highway corresponds to a 
smallest enclosing figure — we will see this theme repeatedly in the following, see Figure [1] Note that the 
result holds in any L p -metric, as all travel to and from the highway is parallel to the y-axis. 
The theorem generalizes to highways of arbitrary orientation in the Euclidean metric: 

Theorem 2 Given n points in the plane, the middle line of the smallest enclosing strip is an optimal 
highway of infinite speed. It can be computed in O(nlogn) time. 

The algorithm used here is the rotating calipers algorithm [13] . After computing the convex hull of the 
point set, it runs in linear time. 



3 The optimal horizontal highway in the Li-metric 

Let p ^ 1 be a real. We say that a p-rhombus is a rhombus of aspect ratio p. 

Theorem 3 Given n points S in the plane, the horizontal axis of symmetry of the smallest enclosing 
axis-aligned v-rhombus is an optimal horizontal highway for the L\-metric and highway speed v. It can 
be computed in 0(n) time. 

Proof. For a pair of points p,q <E S, let d(p, q) := \x p — x q \ / v + \ y p — y q \. Clearly, d(p, q) is a lower bound 
for the travel time diameter of S for any horizontal highway, and therefore 8 := max P:9e s d(p, q) is also 
a lower bound. We show that in fact this bound can be obtained, resulting in an optimal highway. 

We observe that the point set can be enclosed in a rhombus with horizontal diagonal 6v and vertical 
diagonal 5. For an example with v = 2, see Figure^ If we place a horizontal highway along the horizontal 
diagonal of this rhombus, then any point in the rhombus has travel time at most 5/2 to the center of the 
rhombus. This implies that the travel-time diameter is at most S. 

The computation boils down to computing minimum and maximum y-axis intercepts among all lines 
through points of S of slope 1/v and —1/v. ED 



4 The optimal highway in the Euclidean Metric 

In this section we consider the Euclidean metric in the plane, and a highway of finite speed v > 1. 

Theorem 4 Given n points S in the plane, an optimal horizontal speed-v highway under the Euclidean 
metric can be computed in 0(n) time. 
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Fig. 3: Optimal highway under the Z/2-metric. 



Proof. The quickest path (that is, the path with shortest travel-time) between two points p and q has 
one of two forms [2]: It is either the segment pq; or a path consisting of three segments pp + , p + q~ , q~ q, 
where p + and q~ are points on the highway, and the lines pp + and qq~ form an angle of a = arccos 1/v 
with the highway, see Figure I5al 

Now let us define a norm r](x,y) on R 2 as 

rj(x,y) = \x\ cosa + \y\ sina. 

Since < a < tt/2, we have rj{x, y) > unless (x, y) = (0, 0), and rj is indeed a norm. 

Let p and q be two points such that the highway is inbetween p and q and such that the shortest 
path between p and q makes use of the highway. Then the travel time from p to q is rj(q — p). Indeed, 
let (x,y) = q — p, and assume x,y ^ 0. Then the travel time from p to q is ay + (x — by)/v, where 
a = 1/sina and b = 1/tana = a/v. Using 1 — 1/v 2 = 1/a 2 , we get ay + (x — by)/v = n{x,y). 

When the highway cannot be used because the line pq forms an angle larger than a with the highway, 
then the travel time is simply the Euclidean distance d(p, q), and r/iq — p) is an underestimate. 

The unit circle under the norm 77 is a (sin a/cos a)-rhombus, i.e., a (tana)-rhombus. We can find 
the smallest such rhombus enclosing a given set S of n points in linear time. This means we have the 
smallest factor 8 > such that the S fits in the rhombus R with corners (0, 8a), (0,—8a), Sv,0), and 
(— 5v,0) (after translating the point set). 

We claim that the x-axis is now an optimal highway. We already know that there is a pair of points 
whose 77-distance is 26, so this is a lower bound on the diameter. We now show that for any pair of 
points, either their travel time (with respect to the highway at y = 0) is at most 26, or they cannot use 
any horizontal highway. 

For any two points p, q in the rhombus R, we have rj(q — p) ^ 28. This means that if the highway 
lies inbetween, then we are already done. So assume that both p and q lie above the highway, and such 
that they can use a horizontal highway. Let q' be the reflection of q around the highway. Since R is 
symmetric with respect to y = 0, q' is also in R, and if p, q' can use a horizontal highway, then their 
travel-time distance is at most 28, implying that the travel time from p to q is also at most 28. 

It remains to consider the case that p, q' cannot use a horizontal highway. This means that the line 
pq 1 forms an angle larger than a with the x-axis, see Figure [3bl Note that y(q' — p) still has a geometric 
meaning: There is a path from p to the highway, then backwards along the highway, then straight to q, 
see Figure I5cl 

The 77-distance measures the whole travel time, but counting the time on the highway negative. Re- 
flecting the last segment of this path back around the highway, we obtain a path from p to q with 
travel time 77(9' — p), still counting time spent on the highway negative. But now observe that this 
path self-intersects in a point x, see Figure I3dl Let x' be the reflection of x. Then r)(q' — p) = 
d{p, x) + n(x' — x) + d{x' , q') = d{p, x) + T)(x' — x) + d(x, q) ^ d{p, x) + d(x, q) ^ d{p, q) (using i] ^ 0). It 
follows dip, q) < 77(9' - p) < 28. H 
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Fig. 4: Unit disk in the travel-time metric (shaded) and smallest enclosing (tan a)-rhombus. 

Again we found an optimal highway by computing a minimal enclosing shape. Interestingly, the shape to 
be minimized is not the unit circle under the travel-time metric. If the highway is the x-axis then the unit 
circle in the travel-time metric is the convex hull of the points (0, v) and (0, —v) and the Euclidean unit 
circle centered at the origin, see the shaded region in Fig. 3J The reader may wish to verify that finding 
the smallest copy of this (unit-disk) shape enclosing S does not give the optimal horizontal highway. The 
rhombus we are minimizing instead is the "rhombus approximation" of this shape, see Fig. 0] 

The very simple linear-time algorithm above results in a horizontal highway that minimizes the 
travel-time diameter of the point set, but it does not actually tell us what the travel-time diameter is. 
Surprisingly, it is actually not possible to compute the travel-time diameter within the same time bound, 
as the following lemma shows. 

Lemma 1 In the algebraic decision-tree model, computing the travel-time diameter for a set of n points 
and a given highway takes £l(nlogn) time. 

Proof. The following problem has a lower bound of Q(n\ogn) in the algebraic decision-tree model: 
Given two sets A and B of n real numbers, is A n B = 0? We show how to transform this problem in 
linear time into a decision instance of the diameter problem. 

Our instance consists of a set A' of n points and a set B' of n points, computed from A and B. 
All points lie on the unit circle. We first scale all numbers in A and B so that they are close to zero 
(depending on v). For each a in A, we create the point (a, vl — a 2 ). For each b in B, we create the point 
(—6, — Vl — b 2 ). Note that since the points are close to the y-axis, no horizontal highway can be used 
to speed up the connection between A' and B' , and so the diameter of the set is simply the Euclidean 
diameter. It follows that the diameter of A' LIB' is 2 if and only if A and B contain a common number. ED 

Theorem 5 Given a set S of n points in the plane, the optimal highway with speed v can be found in 
0(n 2 logn) time. 

Proof. First we compute the convex hull C of the point set S. Then we use the rotating-calipcrs 
algorithm [12] to compute the function w : [0, w) — > R + that maps an angle <j> to the width of the 
smallest strip that contains C and makes angle cf> with the positive x-axis. The function w consists of at 
most n pieces, each of which is a trigonometric function that can be computed explicitly in constant time. 
Let 7 = 7r/2 — a be the angle formed by the main diagonal and the sides of a (tan a)-rhombus (where 
cos a = 1/v), see Fig. Then the function u((f>) = max{u>((/) — 7 mod 7r), w{4> + 7 mod 7r)} maps (f> to 
the width of a smallest (tan a)-rhombus that contains C and whose main diagonal forms an angle of <p 
with the positive a;-axis. 

As we saw above in the proof of Theorem 2J the main diagonal of this rhombus is an optimal highway 
for this orientation, and the travel-time diameter z{<p) for this highway is 

z{4>) = max{u(</>), max d(p, <?)}, 

where the maximum is taken over all pairs p,q G S that cannot use any highway with orientation <fi. 
This is the case if the vector p — q makes an angle larger than a with the highway orientation. We define 
now, for every pair p, q £ S, a function f pq : [0, n) — > WLq with f pq (4>) = d(p 7 q) if the angle between the 
vector p — q and the positive a;-axis is inbetween (f> + tt/2 — a and <f> + 7r/2 + a modulo tt, and f pq {4>) = 
otherwise. We can then rewrite the optimal travel-time diameter for orientation <f> as 

z{<j>) = max{M(0), max f pq (4>)}. 

p,q€S 
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The graph of z is the upper envelope of (™) horizontal segments and of the graph of u (which has 0{n) 
breakpoints). Thus the graph of z can be computed in time 0(n 2 log n), and we find an optimal orien- 
tation by picking a lowest point on this graph. ED 



The function u used in the proof here can be computed in O(nlogn) time, and within the same time 
bound we can pick an orientation cj) for which u is minimal. This means that u(<fi) is a lower bound for 
the travel-time diameter of any highway for S. Let R be the smallest enclosing (tana)-rhombus with 
orientation cj). Then any point in R has travel-time distance at most u(4>)/(2 sin a) from the center of 
R, see Fig. |U This implies that the highway /i app through the main diagonal of R has a travel-time 
diameter of at most u((f>)/ sin a, which is within a factor of 1/sina = \J v 2 /(v 2 — 1) from optimal. On 
the other hand, any highway yields at least the same (factor-z;) approximation as building no highway 
at all. It is easy to see that the maximum of the function min{w, yjv 2 j(v 2 — 1)} for v ^ 1 is attained at 
v = \/2 and has a value of \/2. Thus ft. a pp is in fact at least a \/2-approximation of the optimal highway. 

Theorem 6 Given n points S in the plane, a speed-v highway with travel-time diameter within a factor 
o/min{t>, \/v 2 1 (v 2 — 1)} ^5 \/2 from optimal can be found in 0(n\ogn) time. 

The approximation factor of this algorithm depends on v, and since 



the factor tends to 1 with growing v very quickly. For instance, for speeds v = 2, 3, and 10, the factors 
are at most 1.16, 1.06, and 1.005, respectively. 

5 The optimal highway cross for infinite speed 

Now we consider the problem of placing more than one highway. Observe that multiple parallel highways 
with the same speed do not reduce the maximum travel time because the quickest path using several 
highways can be simulated with only one highway. Instead we investigate highway crosses, that is, pairs 
of highways that intersect perpendicularly We give algorithms for computing the optimal axis-aligned 
highway cross. 

Definition 1 An enclosing cross for a point set S is the union of a horizontal and a vertical strip of 
equal width containing S . 

Lemma 2 The travel-time diameter of the optimal axis-aligned highway cross with infinite speed equals 
the width of the smallest enclosing cross. 

Proof. Let 6 be the travel-time diameter, and let 8' be the width of a smallest enclosing cross C. 

We first show that 8' ^ 8: Let hi, ^2 be a pair of optimal highways. We assign each point in S to its 
closest highway so that S is partitioned into two subsets: one consisting of points closer to the horizontal 
highway and the other consisting of points closer to vertical highway. We put around each highway the 
narrowest strip containing all the points assigned to the highway. Then both strips have width at most 
(5, otherwise there are two points in the wider strip whose travel-time distance is larger than 8. Therefore 
we can obtain an enclosing cross of width 8 by widening each strip until its width becomes 8. Since 5' 
was minimal, we have 8' ^ 8. 

It remains to show 8 ^ 8': We place highways in the middle of each strip of C. This results in a pair 
of highways with travel-time diameter at most 5' . Since 8 is optimal, we have 8 ^ 5'. ED 

Note that once again the optimal facility corresponds to a minimal enclosing shape. This shape can 
be computed efficiently. 

Theorem 7 Given n points in the plane, the optimal axis-aligned highway cross for infinite speed cor- 
responds to the smallest enclosing strip cross. It can be computed in 0(n log n) time. 

Proof. The characterization follows from Lemma [5J The smallest enclosing cross of a set S of n points 
can be found as follows. 




() 



1. We presort the points by their x- and by their y-coordinatcs. 

2. For a given width oj > 0, we can decide in linear time whether an enclosing cross of width u> exists. 
If it is the case, the enclosing cross can be found in the same time. Our decision algorithm is as 
follows. We slide a vertical strip V of width u> across the point set from left to right. We maintain 
a horizontal strip H of smallest width containing all the points not in V. For each point entering 
V from the right or leaving V from the left, we update H accordingly. If the width of H ever 
becomes uj or less, we answer "yes" and report an enclosing cross. Otherwise, we answer "no" . 

3. The width of the smallest enclosing cross is in the list of numbers L = L x U L y , where L x — 
{xj -Xi\l^i<j^n} and Xi ^ x-x ^ ■ ■ ■ ^ x n is the sorted sequence of x-coordinates. The list 
L y is defined analogously based on the sorted sequence of y-coordinates. 

4. Consider the matrix A with -A[z,j'] = Xj — x n —i+i. The rows and columns of A are sorted in 
ascending order. Using the technique of [10] . we can determine the fc-th element of such a sorted 
matrix in 0(n) time without constructing A explicitly. This gives us a way to do binary search on 
L x . This search consists of 0(log n) steps, each of which first invokes the algorithm of Fredcrickson 
and Johnson to find the median of the remaining elements in L and then calls the decision algorithm 
of stage [5J Thus, the total runtime is 0(n log n). Likewise we can search for the smallest value in 
L y for which an enclosing cross exists. Finally we return the minimum of the two values. Notice 
that the algorithm in stage [2] computes not only the width of the smallest enclosing cross but also 
the cross itself. ED 

6 The optimal axis-aligned highway cross for finite speed 

In this problem we have been unable to characterize the optimal solution by a smallest enclosing shape. 
There is no natural "center" for the problem: sometimes there are critical paths where a point connects 
to the highway that is further away. 

Theorem 8 The optimal axis-aligned speed-v highway cross under the L\-metric can be computed in 
0(n a(n)) time. 

Proof. The optimal solution corresponds to the lowest point on the upper envelope of the pairwise 
distance functions, see Figure [5] Since these functions are piecewise linear and of constant complexity, 
their upper envelope can be computed in 0(n 4 a(n)) time [5]. EQ 



6.1 The decision problem 

We now present an algorithm that decides for a given 5 > whether there is an axis-aligned speed-w 
highway cross such that the resulting travel-time diameter is at most S. This will be used as a subroutine 
for finding approximations of the optimal highway cross, see Section [7] 

Theorem 9 Given a set P of m pairs of points in the plane, a speed v > 1, and a parameter S > 0, we 
can decide in time 0(ma(m) log m) whether there is an axis-aligned highway cross of speed v such that 
the travel-time distance of all pairs in P is at most S. If the answer is positive, we can compute such a 
highway cross within the same time bound. 

Proof. For points o~,p,q 6 M 2 , let d a (j>, q) denote the travel-time distance between p and q, assuming 
an axis-aligned highway cross with speed v has been placed (with center) at a. We define the region 

R(p,q) :={a eM. 2 \d a {p 7 q) ^S}. 

We observe that the answer to the decision problem is positive if and only if ("L g)gp R{Pi o) i s n °t 
empty 

The shape of the region R(p,q) depends on 5. Let w,h be the horizontal and vertical distance of 
points p and q. If 8 < (w + h)/v, then R(p, q) is empty If (w + h)/v ^ 6 < min{w + h/v,h + w/v}, then 
R(p, q) consists of two convex quadrilaterals. If min{u; + h/v,h + w/v} ^ 6 < max{w + h/v, h + w/v}, 
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Fig. 5: Travel-time distance between two points p and q for an axis-aligned highway cross centered at {x,y). 

The graph of this function has 31 faces of 15 different orientations. The two points marked x are the 
lowest, that is, those where the corresponding highway crosses minimize the travel time from p to q. The 
thin dark gray lines are contour lines. If a highway cross is centered in a V- or H-region, the vertical 
and horizontal highway, respectively, is used by a quickest p—q path. In the VH- and HV-region both 
highways are used in the corresponding order. In the N-regions no highway is used. 







F A (p,q) F 3 (p,q) 



(P,q) F 2 (p,q) 




(a) 



(b) 



Fig. 6: (a) the regions R(p, q) when min{™ + h/v, h + w/v} ^ S < max{ui + h/v, h + w/v} (light gray region) 
and when max{w + h/v, h + w/v} ^ 6 < w + h (dark and light gray regions), (b) the dark and light 
gray regions can be expressed as the intersection of the four types of regions F\,F2,F3, and F4. 
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Fig. 7: Here the optimum speed-u highway cross (dashed) is not contained in the optimum speed-oo strip cross 



then R(p,q) is infinite in one (axis-parallel) direction. If max{w + h/v,h + w/v} ^ S < w + h, then 
R(p, q) is infinite in both axis-parallel directions. Finally, if w + h ^ 6, then R(p, q) = R 2 . See Figure [6l 
Let us call a planar region F (a, 6)-monotone if for every point {x, y) G F and any A ^ the point 
(x + Xa, y + Xb) is also in F. We observe that R(p, q) can be expressed as the intersection of four re- 
gions Fi(p,q), i = 1,2,3,4, where F\{p,q) is (1, l)-monotone, F^{p, q) is (— 1, l)-monotone, F 3 (p 7 q) is 
(— 1, — l)-monotone, and F±(p, q) is (1, — l)-monotone. Figure El shows an example of R(p, q), which can 
be expressed as the intersection of Fi,F 2 ,F 3 , and F4. Each region is bounded by a polygonal curve of 
constant complexity. F, := |"), , £ p Fi(p, q) is the lower envelope of a set of 0(m) line segments, has 
complexity 0(ma(m)) [T2], and can be computed in time 0(m\ogm). The intersections F\ n F3 and 
F 2 n F4 can be computed by a plane sweep in time 0{ma{m) log m). We are left with two regions of 
complexity 0(ma(m)), and we need to determine whether their intersection is empty. While we do not 
know how to bound the complexity of this region, we can test emptiness in O (ma (to) log to) time, by 
a simple plane sweep that stops as soon as a point in the intersection is found. Since any intersection 
between edges of the two regions implies that the intersection is not empty, this runs in the claimed time 
bound. If such an intersection is found, it is a center for highway cross with travel-time diameter at 
most a. ED 



6.2 Further observations 

Suppose we could characterize the travel-time diameter given the optimal highway cross to get a compact 
list L of candidate values as in the case of the infinite-speed highway cross. Then we could do binary 
search on L using the decision algorithm of Theorem [5] 

Given the travel-time graph T pq for each pair of points p and q (see Figure [SJ , consider the upper 
envelope over all these graphs as in Theorem [8] The minimum occurs at a vertex of this upper envelope. 
Thus there are always at most three pairs — that is, at most six points — that define this vertex. Locally, 
in a neighborhood of the optimal highway, these at most six points alone have the same optimal highway 
cross. However, the lower envelope for these three pairs might have other minima, and the original 
solution might not be the global minimum. We do not know whether this is the case, and whether it 
could perhaps be resolved by considering a few more pairs. If that was possible, than we could apply 
Chan's technique [8] (as Cardinal and Langerman do [7]), in order to get a randomized algorithm for 
the optimization problem whose expected running time is asymptotically the same as the running time 
of the decision algorithm. 

The optimum axis-aligned highway cross for finite speed need not be contained in the strip cross 
for infinite speed, see Figure [7] take the points (-2,1), (-1,2), (1,2), (2,1), (2,-1), (1,-2), (-1,-2), 
(—2,-1) — that is, an octagon, contained in the strip cross ([—1,1] x 1) U (1 x [—1,1])) — plus the two 
points (—a, 0) and (a, 1) for a > 3w. Then the coordinate axes are the optimum highway cross for infinite 
speed. It has travel-time diameter 2, so the above two strips are the optimum cover. But for any finite 
speed v > 1, the highway cross x = a and y = yields a diameter of (2a + l)/v, which is better than 
the diameter (2a/v) + 1 caused by the coordinate axes being highways. 

The following argument also rules out any simple incremental algorithm: Even for infinite speed there 
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Fig. 8: After adding the point (3, 10) the optimal speed-oo highway cross changes (from solid to dashed), but 
the new point does not occur in any diametral pair. 

are point sets such that the addition of one point changes the diameter, and the new point does not 
occur in any diametral pair. An example (see Figure [5]) for infinite speed is given by the points (0,9), 
(0, 11), (9, 0), (11, 0) since now the coordinate axes are an optimal highway cross — with a diameter of 0. 
If we add the point (3, 10), the optimal highway cross is centered at (10, 10) and has a diameter of 2. 

7 Approximations for the optimal axis-aligned highway cross 

Given a set S of n points. Let C be the smallest enclosing cross for S, and let hi, h^ be the middle line 
of each strip of C . We call hi, /12 the median highways for S. 

Lemma 3 The travel-time diameter S mc d of the median highways ( with speed v) is at most 2 + 1/v times 
the travel-time diameter 6 op t of an optimal axis-aligned speed-v highway cross for S . There are point sets 
S where for v \/3 the travel-time diameter of the median highways is at least 2 — l/(v + 2) times the 
optimum. 

Proof. We can scale S such that its Li-diamctcr is 2 — this does not change the travel-time ratio. Let 
w be the width of C after scaling. Observe that 5 op t ^ 2/v, as there are points at Li-distancc 2. 
Furthermore, we have <5 op t w, since using the optimal highways at infinite speed cannot achieve 
diameter less than w. 

On the other hand, <5 mo d w + (2 + w)/v, since any point can reach a point on the highways at 
distance at most w/2, and the maximum distance of such points on the highways is at most 2 + w. This 
implies <5 mcd < (1 + l/v)w + 2/v < (1 + l/v)S opt + S opt = (2 + l/v)S opt . 

For the lower bound example, let speed v > 1 be given, and set parameter u> = l/(v + 2). We will 
construct a point set S such that the smallest enclosing cross has width 2lu, the median highway has 
travel-time diameter 4/v — 2/(v(v + 2)), and the optimal highway cross has travel-time diameter at most 
2/v, implying the lower bound. 

Let e > be very small. Our point set S consists of the points (1,0), (—1,0), (0,1), (0,-1), 
(— 2oj, uj + e), (— 2uj, — uj — e), (e, — 2u>), as in Figure [9l We claim that S has a unique smallest enclosing 
strip of width 2uj, centered around the lines x = —uj and y = —u>. Indeed, the line x — must be in the 
vertical strip (otherwise the horizontal strip would have width at least 2), while the line y = must be 
in the horizontal strip. Similarly, the line x = —2uj must be in the vertical strip as well, and this now 
fixes the vertical strip of width 2u> around the line x = —uj. It follows that the remaining point (e, — 2u>) 
is in the horizontal strip, fixing that strip around y = —uj. 

The median highway cross has travel-time diameter 2uj-\-{2+2uj) / v = 4/v— 2/(v(v+2)) if the diameter 
is less than or equal to the Li distance of these points (note that the diameter is determined by (1,0) 
and (0, 1).) That is, for v ^ -\/3 the median highway cross has travel-time diameter A/v — 2/(v(v + 2)). 
Consider now a highway cross with center at the origin. The four outer points and the point (e, —2uj) 
can be reached from the origin within travel time 1/v. The remaining two points can be reached with 
travel time uj(1 + 2/v) (ignoring all e-terms). The travel-time distance between these two points is 2uj, 
and so the travel-time diameter is bounded by 



max{2/w,l/w + u(l + 2/v),2uj} = 2/v. 



□ 
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Fig. 9: The median highway for u = 2. 

We can improve the result in Lemma [3] by a simple observation. 

Theorem 10 Given a set S of n points we can compute in 0(n log n) time an axis-aligned highway cross 
whose travel-time diameter is at most 1 + \[2 times the travel-time diameter of an optimal axis-aligned 
speed-v highway cross for S. 

Proof. According to Theorem [7] the median highways can be computed in 0(?ilogn) time. According 
to Lemma [3] they yield a factor-(2 + 1/v) approximation for the optimal travel-time diameter. Note that 
the approximation factor tends to 3 when the speed goes to 1 . Clearly not building a highway cross is a 
factor-w approximation. Balancing out the two terms yields min{2 + 1/v, v} ^ 1 + \/2. ED 

The following lemma will allow us to do better. 

Lemma 4 Let s be any point in S . Let H s be the highway cross that minimizes the maximum travel 
time to s. The travel-time diameter of H s is at most twice the travel-time diameter S op t of an optimal 
axis-aligned speed-v highway cross for S. 

Proof. Let {p, q} be a pair of points in S and let s' G S be a point of maximum travel-time distance 
from s given H s . We denote by d s the metric induced by H s and by d op t the metric induced by the 
optimal highway cross. Then d s (s,p) ^ d s (s,s') ^ d op t(s,s') ^ <5 pt and, by symmetry, d s (s,q) ^ <5 opt . 
This yields d s (p, q) sC d s (p, s) + d s (s, q) ^ 2<S opt . ED 

Note that H s is usually not centered at s (consider, for instance, the set S = {(0, 1), (1,0)} whose 
optimal highway cross is centered at the origin). 

Based on the constant-factor approximation from Theorem \W\ we can use binary search and the 
decision procedure of Theorem H] to get the following. 

Theorem 11 Given a set S of n points in the plane, we can compute in 0(log(l/e)a(n)nlogn) time a 
(2 + e) -approximation for the optimal axis-aligned speed-v highway cross for S . 

Proof. Let s be any point in S and H s be the highway cross that minimizes the maximum travel time 
to s. 

According to Theorem|5]the travel-time diameter <5 S given H s is at most twice the travel-time diameter 
<5opt given an optimal axis-aligned spced-u highway cross for S, that is. 5 S ^ 2J opt , so a (1 + e/2)- 
approximation 5 £ to S s is a (2 + ^-approximation for S opt - 

We now describe how to compute a (1 + e)-approximation for S s by binary search. Recall that the 
median highways yield a travel-time diameter of <5 mc d ^ (l + V / 2)<5o P t ^5 3<5 op t ^ 35 s oiS s , see Theorem [TU1 
The median highways can be computed in O(nlogn) time according to Theorem [7] 
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Now we conceptually subdivide the interval I = [0, 2(5 me( j] into at most TV = 6/e pieces of length 
5 mcc i • e/3, and denote the increasing sequence of interval endpoints by A = (Si, . . . ,Sn). Since S s ^ 
2i5 pt ^ 2(5 me d, we know that S s lies in /. Hence there is an index i G {1, . . . , TV} such that Si < S s ^ tfj+i. 

Setting S e = Si+i we find that S s ^ Si = S e —S me d-s/3. This yields S £ <5 s + <5 mc d-e/3 S s + 3S s -e/3 ^ 
(1 + e)S s . Thus S s is indeed a (1 + e)-approximation of S s . 

For a given 8 > we run the decision algorithm of Theorem [5J using the set of n — 1 pairs 
P = {(s,q) \ q E S \ {s}}. Each such test takes 0(a(n)nlogn) time. Using 0(log(l/e)) calls to 
this decision procedure, we can determine S e by binary search on A. We return the highway cross com- 
puted by the decision procedure for the largest Si ^ S s . ED 

If we are willing to invest more time, we can even get a (l+e)-approximation of the optimal travel-time 
diameter S op t- 

Theorem 12 Given a set S of n points in the plane, we can compute in 0(log(l/e)a(n)n 2 \ogn) time 
a (1 + e)- approximation for the travel-time diameter of S under the optimal axis-aligned speed-v highway 
cross. 

Proof. We again first compute the median highways to get an upper bound S me d for the optimal travel- 
time diameter <5 opt and then do binary search. We can now use the interval [0,<5 mc d], which contains 
Sopt- We stop when the interval size is sufficiently small, that is, at most i5 mo d • e/3. This time we use 
the decision algorithm of Theorem [5] with the set P of all (™) pairs of points in S. ED 



8 Concluding remarks 

There arc many ways how this problem can be extended. First, can we compute an optimal highway with 
arbitrary orientation under the Euclidean metric in o(n 2 ) (worst-case) time? Second, consider highways 
with different speeds, different slopes, or bounded lengths. Third, suppose an existing network of (axis- 
parallel) highways and a real I > is given. Where to place a new (axis-parallel) highway segment of 
length I in order to minimize the travel-time diameter of the resulting network? 
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